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Abstract 

We study a two-dimensional family of probability measures on infinite 
Gelfand-Tsetlin schemes induced by a distinguished family of extreme 
characters of the infinite-dimensional unitary group. These measures are 
unitary group analogs of the well-known Plancherel measures for symmet- 
ric groups. 

We show that any measure from our family defines a determinantal 
point process on Z+ x Z, and we prove that in appropriate scaling limits, 
such processes converge to two different extensions of the discrete sine 
process as well as to the extended Airy and Pearcey processes. 



1 Introduction 

Let <S'(7i) be the symmetric group of degree n. Denote by Y„ the set of partitions 
of n or, equivalently, the set of Young diagrams with n boxes. It is well known 
that complex irreducible representations of S{n) are parameterized by elements 
of Y„; we denote by dim A the dimension of the irreducible representation cor- 
responding to A. The probability distribution 

Prob{A}=^^^, AeY„, 

on Y„ is called the Plancherel measure for S{n). The Plancherel weight of 
A G Y„ is the relative dimension of the isotypic component of the regular repre- 
sentation of S{n), which transforms according to the irreducible representation 
corresponding to A. Hence, one has the following equality of functions on S{n): 

^ _ dim^ A x'^ 
^ ^-^ n\ dim A 

AG¥„ 

where (5e is the delta-function at the unity, and x^ is the irreducible character 
corresponding to A. 

Let S'(oo) = U„>iS'(n) be the group of finite permutations of a countable set 
known as the infinite symmetric group, see e.g. [17'. The group S{oo) has a rich 



theory of characters (positive-definite central functions on the group). For any 
character x of <S'(oo) normalized by x(e) = 1, its restriction to the subgroup S{n) 
of permutations of first n symbols is a convex combination of {x^ / dimA}AeY„- 
The coefficients x„(A) form a probability measure on Y„; they are a kind of 
Fourier transform of x- 

There exists only one character % of 5(oo) for which the rows and columns of 
the Young diagrams distributed according to x„ grow sublinearly in n as n ^ oo. 
This character is the delta-function at the unity of 5'(oo), the corresponding 
representation is the (bi)regular representation of ^(oo) in ^^(5(00)), and Xn is 
the Plancherel measure on Y„ introduced above. 

An analogous construction for the infinite-dimensional unitary group U (00) = 
^N>iU{N) yields a two-dimensional family of characters of ?7(oo). Although 
the notion of regular representation for U{oo) is meaningless, by comparing the 
lists of the extreme (i.e., indecomposable) characters of 15(00) and U{oo) one 
sees that the analog of 5^ on 5(oo) is the family of characters 



([/) = exp (Tr (7+(C/ - 1) + 1~{U~^ - 1))) , U G C/(oo), 



where 7^ > are the parameters of the family. We will provide details in 
Section 3, and for now let us just say that on the level of Fourier transform, 
the set Y„ is replaced by the set of A^-tuples of integers Ai > • • • > Aat which 
we call signatures or highest weights of length A'' (they parameterize irreducible 
representations of the unitary group U{N)), and the corresponding probability 
distributions have the form 

P^^'^"(Ai, ...,Xn)= const-det[/f ^'^"^(A, - j)]J',-^idimc;(^)(A), 

where dim[/(jv)(A) is the dimension of the irreducible representation of U{N) 

with highest weight A. We call the measures the Plancherel measures 

for the infinite- dimensional unitary group, and the present paper is devoted to 
the study of these measures. 

One source of interest to the Plancherel measures for symmetric groups is 
the fact that the distribution of the largest part of A € Y„ coincides with 
the distribution of the longest increasing subsequence of uniformly distributed 
permutation in S{n). This fact can be restated in terms of a random growth 
model in one space dimension called the polynuclear growth process (PNG). 
Namely, the distribution of the height function for PNG with the so-called 
droplet initial condition at any given point in space-time coincides with the 
distribution of the largest part of A € U„>oY„ distributed according to the 
Poissonized Plancherel measure 

6(1^1 dim A' ^ 



Prob{A} = e-*' ( ) , AeU„>oY„, 



2 



where |A| is the number of boxes in the Young diagram A, and > is a 
parameter, see [25] , 

Quite similarly, the largest coordinate of a signature distributed according 
to the Plancherel measure for U{oo) describes the height function in another 
growth model in one space dimension called PushASEP for the so-called step ini- 
tial condition. This fact can be established by direct comparison of Proposition 
3.4 from [B] and Theorem 3.2 below. 

The asymptotics of the Plancherel measure for S{n) as n ^ oo has been 
extensively studied. In the seventies, Logan and Shepp [19] and, independently, 
Vershik and Kerov [55], [3D], discovered that Plancherel distributed Young di- 
agrams have a limit shape: In a suitable metric, the measure on these Young 
diagrams scaled by ^/rl converges as rt ^ oo to the delta-measure supported 
on a certain shape. In the late nineties, more refined results were obtained. It 
was shown that the random point process generated by the rows (or columns) 
of the Plancherel distributed Young diagrams has two types of scaling limits, in 
the "bulk" and at the "edge" of the limit shape. In the limit, the former case 
yields the discrete sine determinantal point process, while the latter case yields 
the Airy determinantal point process, see [3], [2], [H], [1], [Hj. 

The main goal of the present paper is to prove similar asymptotics results 
on scaling limits of random point processes related to more complex measures 
with N ^ oo and 7* possibly dependent on A^. Note that our results 
do not imply the existence of the limit shape in any of the cases we consider, 
although they strongly suggest that in some cases the limit shape does exist, 
and they predict what it looks like. For a discussion of the relationship between 
"local" results on point processes and "global" measure concentration properties 
see Remark 1.7 of [S], §1 of [12] . 

Let us describe our results in more detail. 

It is convenient to represent a signature A — {Ai > • • • > Xn} as a pair of 
partitions, one partition A"*" consists of positive parts of A while the other one 
A~ consists of absolute values of negative parts of A. When the parameters 7* 
are independent of A^, they describe (see Section 2) the asymptotic behavior of 
|A^|, namely |A^| ~ 'J^N, as N ^ 00. This asymptotic relation remains true in 
other situations as well, and it is helpful to keep it in mind when going through 
the limit transitions below. 

Our first result describes what happens when 7* ~ A^^ as A ^ 00. Then 
one expects that jA"*"!, |A~| remain finite in the limit, and indeed the measures 

'"^ ■* converge to the product of two independent copies of the Poissonized 
Plancherel measures for the symmetric groups that live on A*. 

The next possibility to consider is when 7* are independent of A^. The 
case when 7^ = was considered by Kerov [TB], who proved the existence 
of the limit shape and showed that the limit shape coincides with that for the 
Plancherel measures for symmetric groups. We show that when both parameters 
7* are fixed and nonzero, the random point processes describing A* asymptot- 
ically behave as though A* represent two independent copies of the Poissonized 
Plancherel measures for the symmetric group with Poissonization parameters 
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7=^iV oo. 

The most interesting case is when 7^ grow at the same rate as N. Biane 
|4] proved that when 7" = 0, the corresponding measure has a hmit shape that 
depends on the Hmiting value of the ratio 7+/iV. We consider the case when 
both parameters are nonzero and investigate the asymptotic behavior of the 
random point process that describes our random signatures. 

Even though we do not prove the existence of the limit shape, it is convenient 
to use the hypothetical limit shape inferred from the limit of the density function 
to describe the results. There are three possibilities: The limit shapes of 
scaled by N do not touch (that happens when 7^/iV are small), when they 
barely meet, and when they have already met, see Figure |4] in the body of the 
paper. Accordingly, there are three types of local behavior one can expect: The 
bulk, the edge, where the limit shape becomes tangent to one of the axes, and 
the point when the edges of the limit shapes for meet. We compute the 
local scaling limits of the correlation functions for the random point process 
describing our signatures, and obtain the correlation functions of the discrete 
sine. Airy, and Pearcey determinantal processes in the three cases above. 

As a matter of fact, we consider probability measures on a more general 
object than signatures. Every character of C/(oo) naturally defines a probabil- 
ity measure on Gelfand-Tsetlin schemes (a kind of infinite semistandard Young 
tableaux) , see Section 2 and references therein. The corresponding measures on 
signatures of length N are certain projections of the measure on Gelfand-Tsetlin 
schemes. In particular, every character from our two-dimensional Plancherel 
family yields a measure on Gelfand-Tsetlin schemes, and that is what we study 
asymptotically. We interpret each scheme as a point configuration in Z x Z-|- , and 
compute the scaling limits of correlation functions of the arising two-dimensional 
random point processes. The results are appropriate (determinantal) time- 
dependent extensions of the limiting processes mentioned above. 

The proofs are based on the techniques of determinantal point processes. 

First, we show that for any extreme character of U{oo), the corresponding 
random point process on Z x Z+ is determinantal, and we compute the correla- 
tion kernel in the form of a double contour integral of a fairly simple integrand. 



This result (Theorem 3.2 1 is similar in spirit to the formula for the correlation 
kernel of the Schur process from [23l , but it does not seem to be in direct rela- 
tionship with it. After that we perform the asymptotic analysis of the contour 
integrals largely following the ideas of [50] , [55] , [H] . 

Acknowledgements. The authors are very grateful to Grigori Olshanski for 
a number of valuable suggestions. The first named author (A. B.) was partially 
supported by the NSF grant DMS-0707163. 



2 Description of the Model 

Let U{N) denote the group of all x unitary matrices. For each A^, U{N) 
is naturally embedded in U{N + 1) as the subgroup fixing the (A^ -I- l)-th basis 
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vector. Equivalently, each U G U{N) can be thought of as an (iV+ 1) x {N + 1) 
matrix by setting Ui^N+i — Un+i.j — for 1 < i, j < N and Un+i,n+i = 1- 
The union W^^^UiN) is denoted U{oo). 

A character of U{oo) is a positive definite function x '■ U{oo) — s- C which is 
constant on conjugacy classes and normalized (x(e) = !)• We further assume 
that X is continuous on each U{N) C U{oo). The set of all characters of [/(oo) 
is convex, and the extreme points of this set are called extreme characters. 

The extreme characters of U{oo) can be parametrized as follows: Let 
denote the product of countably many copies of M. Let ft be the set of all 
(a+,a-,/3+,/3-,(5+,(5~) such that ([25], §1) 

a± = (af > > . . . > 0) e = {13^ > > ■ ■ ■ > 0) £ (5± e M, 



Set 

CSO 

i=l 

Each in this set defines a function x'^ on U{oo) by 

x"(t/) = n 

uGSpcctrum(C/) 



f^i 1 - al{u - 1) 1 - a, (u 1-1) 
As uj ranges over fi, the functions turn out to be all the extreme characters 

of f/(c3o) (ini, ES], ESj). 

Equipping IR°° x x x x E x E with the product topology induces 
a topology on fl. For any fixed U G C/(oo), x"(t^) is a continuous function of u). 
For any character x of [/(oo), there exists a unique Borel probability measure 
P on ri such that 

x{u)= [ x'^mdp, 

Jn 

see |25j . Theorem 9.1. This measure is called the spectral measure of x- 

It is a classical result that the irreducible representations of U{N) can be 
paramterized by nonincreasing sequences A = (Ai > ... > Ajv) of iV integers 
(see e.g. [32J). Such sequences are called signatures (or highest weights) of length 
N. Thus there is a natural bijection A ^ x^ between signatures of length N and 
the conventional irreducible characters of U{N). 
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The extreme characters of U{oo) can be approximated by with growing 
signatures A. To state this precisely we need more notation. 

Represent a signature A as a pair of Young diagrams (A+, A^), where A+ 
consists of positive A^'s and A~ consists of negative A^'s. Zeroes can go in either 
of the two: 



A — (A^, AJ, . . . , — A2 , — Aj ). 

Let d(-) denote the number of diagonal boxes of a Young diagram and set 
(i+ = d{X'^) and d~ = d{X~). Recall that the Frobenius coordinates Pi,qi of a 
Young diagram A are defined by 

Pi = Ai - z, = (A')j - j, 1 < i < rf(A), 

where A' is the transposed diagram. 

The dimension of the irreducible representation of U (N) indexed by a sig- 
nature A = (Ai, . . . ,Xn) is given by Weyl's formula: 



dim^A = x\l,---,l) = n 



Aj - i - Xj + j 



1 — I 



Define the normalized irreducible characters by 

~^ dimArA'*'" 

Note that x^{e) = 1. 

Given a sequence {/at} of functions on U{N), we say that /at's approximate 
a function / on U{oo) if for any fixed Nq, the restrictions of the functions 
(for > A^o) to U{No) uniformly tend, as — > 00, to the restriction of / to 
U{No). We have the following approximation theorem: 

Theorem 2.1. Let x be the extreme character corresponding to (a*,/3*,7^) € 
f2. Let {A(n)} he a sequence oj signatures oj length n with Frobenius coordinates 



p. 



{n),q^{n). Then the functions approximate x iff 

lim^..?, lim^./3^ liml(^("))^ 



-^00 Tl TL — >oo TL n — >oo 



for all i. 

Proof. This theorem is due to Vershik and Kerov [S3]. See [32] for a detailed 
proof. □ 

Let GTat be the set of all signatures of length N and set GT = UatGTjv- 
Turn GT into a graph by drawing an edge between signatures A G GTat and 
/i S GT^v+i if A and /i satisfy the branching relation A ^ /i, where A ^ 
means that /ii < Ai < /l(2 < A2 < . . . < Aa? < hn+i- GT is also known as the 
Gelfand-Tsetlin graph. 
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Each character of U{oo) defines a probability measure Pn on GTat. If we 
restrict the extreme character to U{N), we can write 

X"k(iv)- E PnWx^- (2) 

AeGTjv 

Definition 2.2. The measure Pn corresponding to the extreme character with 
= (3^ = and arbitrary 7* > will be called the Nth level Plancherel 

measure with parameters 7* . Denote it by P^ . 

The choice of the term is explained by the analogy with the infinite sym- 
metric group S{oo). The extreme characters of S'(oo) are parameterized by 

{(a,/3,7) eM^ xM- xM+;E(". + A)+7 = l}- 

The measure on partitions of n obtained from the character with ai = Pi = 
0,7 = 1, simiarly to the measure Pn above, assigns the weight (dimA)^/n! to a 
partition A and is commonly called the Plancherel measure. Here dim A is the 
dimension of the irreducible representation of Sn corresponding to A. 

Let X be a character of U{oo) and let P and Pn be its corresponding decom- 
posing measures on and GT^r. For any N, embed GTn into CI by sending A 
to (a"*", a^, 6^, c+, c^) where 

± _pf u± - it „± _ l^^^l 



N' ' TV ' N 
Define a probability measure Pjy on f2 to be the pushforward of Pn under 
this embedding. Then Ppf weakly converges to P as TV ^ 00 ([25], Theorem 
10.2). 

This implies that as ^ 00, the Plancherel measures P^ converge to 
the delta measure dX ui — (a^ — (if = 0,7"'", 7"), that is, the row and column 

lengths for A* distributed according to P"^ grow sublinearly in N. 

The main goal of this paper is to study the asymptotic behavior of the 
signatures distributed according to the Plancherel measures P^ as TV ^ 
00. We will also study a more general object: the corresponding probability 
measures on objects called paths in GT. 

A path in GT is an infinite sequence t — (ti, ^2, • • •) such that ti G GT.; and 
ti -< ii+i. Let T be the set of all paths. 

We also have finite paths, which are sequences t — (ti, T2, . . . , tjv) such that 
Ti G GTj and ti ^ T2 ^ . . . ^ tn ■ The set of all paths of length N is denoted 
by Tjv. For each finite path t G T^r, let C,- be the cylinder set 

Cr^{teT: {ti,t2,...,tN)^T}. 

A character x of U{oo) also defines a probability measure on T which 
can be specified by setting 
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where Pn is as above and r is an arbitrary finite path ending at A ([21], §10). 
In particular, any lu € defines a measure on T via the corresponding extreme 
character x'^ ■ If satisfies — (3^^ — with arbitrary 7* , then let P'^ '''' 
denote this measure. 



3 Plancherel measures as determinantal point 
processes 

In order to analyze P"^ and P'^* , it is convenient to represent signatures 
as finite point configurations (subsets) in one-dimensional lattice. Assign to 
each signature A G GT^v a point configuration C{\) C Z by 

A = (Ai, . . . , Ajv) ^C{\) = {Ai - 1, . . . , Aat - TV}. 

The pushforward of under this map is a measure on subsets of Z, 

that is, a random point process on Z. Denote this point process by P]^ ''^ . 
The map A ^ can be seen visually. For example, if A = (4, 2, 0, 0, —1, —3), 
then C{\) = {3, 0, -3, -4, -6, -9}. Sec Figure [l] 



Figure 1: Black dots represent points in the configuration and white dots rep- 
resent points not in the configuration. 
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Given a point process on Z, define the nth correlation function p„ by 

p„ : Z" ^ [0, 1] 
{xi,X2, . ..,Xn) ^ Prob({X C Z : {xi,X2, ■ ■ ■ ,a;„} C X}). 
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(There is a more general definition of correlation functions, but it will not be 
needed here. See e.g. [5], §5 for more details). Clearly, this function is symmetric 
with respect to the permutations of the arguments. 

On a countable discrete state space (Z in our case) a point process is uniquely 
determined by its correlation functions (see e.g. |5], §4), so to study the measure 
it suffices to study its correlation functions. 

A point process is determinantal if there exists a function K such that 

Pnixi,X2, ... ,x„) ^ det[K{xi, Xj)]i<i^j<n for any n = 1,2, 

The function K is the correlation kernel. A useful observation is that K is not 
unique: K{x,y) and j|^i4r(a;, y) define the same correlation functions for an 
arbitrary function /. 

Just as A 1-^ 'C(A) defines a map from GTji to the set of subsets of Z, we 
have a map from the set T of paths in the Gelfant-Tsetlin graph to subsets of 
Z_|_ X Z. Let t = {ti < t2 < . . .) he a, path in GT. Each ti is a signature of 
length i which will be written as A« = (Af\A^^\...,Af ). Then map t to 

C{t) = {(i, A^.'^ - j) : 1 < i < oo, 1 < j < i} C Z+ X Z. 

The pushforward of P''^ ''^ under this map will be denoted by 7'''^'''' . This 
is a random point process on Z+ x Z. 

One more introductory concept is needed. Define a map A by 



Given a point process V on Z+ x Z, its pushforward under A is also a point 
process on Z+ x Z, which will be denoted Pa- The map A is often referred 
to as "particle- hole involution". With this notation, we have the following 
proposition: 

Proposition 3.1. IfV is a determinantal point process with correlation kernel 
K(ni,Xi;nj,Xj), then Va is also a determinantal point prcess. Its correlation 
kernel is ^rLi—nj^xi—xj I^i^rii^ x^^ nj ^ Xj^ . 

Proof. See Proposition A. 8 of ^8]. □ 
Let us now state the main theorem of this section. 

Theorem 3.2. The point process V^^ is determinantal. Let K{ni,Xi;nj,Xj) 
denote its correlation kernel. If ni > n2, then 



^ . 2Z+XZ 



X 



X >^{Z+x Z)\X. 




If ni < 712, then 
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In these expressions, u is integrated over |m| = r < 1 and w is integrated over 
|w— l| = e<l — r and z is integrated over |z| = r < 1. 

Corollary 3.3. The point process V\ is determinantal. LetK/^{ni 
denote its correlation kernel. If ni > 112, then 



KA{ni,xi;n2,X2) 
If ni < n2, then 



27riy 7 7 e'y «'+7+i" 'tyi+^a (1 - w)"^ ^ _ ' 



KaItit ,XT;n2,X2) = <i> -, ^ dz 

^ ' ' ' ' 2-Ki 7 (1 - 

^Yll eT""+T^""'u="i(l dudw 

In these expressions, u is integrated over |w| = r < 1 and w is integrated over 
|w — l| = e<l — r and z is integrated over |z| = ?' < 1. 

Proof. The corollary follows immediately from Proposition |3. 1 1 and the fact that 
^--^^^2;^^,|^, (1-z)— 1^"^ 



for ni ~ 7i2. Note that in Theorem 
and n\ < n2, while in Corollary |3. 3 



3.2| the two cases for the kernel are ni > n2 
the two cases are ni > 712 and ni < ?i2. □ 

Remark. Let K{ni,xi;n2,X2) and KA{ni^xi;n2,X2) denote the correla- 
tion kernels of '^"^ and 7^2 '''^ i respectively (7+ and 7" switched places). 
The substitutions u i— > u~^,w 1— > and further deformation of the contours 
show that 

(-l)"i""^X(ni, -xi - ni - I]n2,-X2 - n2 - 1) = K{ni,xi\n2,X2)^ 
(_l)"i-"2^^(-^^^ -xi - ni - I;n2,-X2 - 71.2 - 1) = KA{ni,Xi; n2, X2). 

This can be understood independently. Switching 7+ and 7^ corresponds to 
switching A"^ and in a signature A. In terms of C{X), this corresponds to 
replacing Xi with —Xi — — 1. For example, consider A = (4,2,0,0,-1,-3) 
from Figure [1] Switching A+ and A" gives A = (3, 1, 0, 0, -2, -4). Then £(A) = 
{2, —1, —3, —4, —7, —10}, which can be obtained from £(A) by replacing Xi with 
—Xi — 6 — 1. 

Remark. The arguments below actually prove a more general statement. 
If we define a point process of Z+ x Z similarly to V •'^ , but starting from an 
extreme character of U{oo) with arbitrary parameters (a^, /3j^, 7^), then this 
process is determinantal and its kernel has a similar form. The only change is 
replacing E{z) below by fo{z) from equation ([ij. 

In what follows we use the notation 

£;(z) = e7+(.-i)+7-(--^-i) = 
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Lemma 3.4. Suppose X = (Ai, A2, . . . , Xn) G GTjv- Write xt for Xk — k. Then 
Pn (A) = const • det[/j(xfe)]i<j,fe<jvdet[5rj(xfe)]i<j,fc<Ar 

where 



fj{xk) = ^ </" E{u)u-^-'''-'du, l<j<N. (6) 

g,{xk)=xi-\ 1<3<N. (7) 

00 

Proof. Writing E{u) ~ c(Z)u' and integrating E{u)u^ over the unit circle, 
/— — 00 

we can solve for c to get 

c{l) = ^ I E{u)u~^^Uu 

Set Lu = {af = (3f = 0,7+, 7"). For U € U{N) with spectrum ui, . . . ,un, 
we can write x'^(^) = E{ui) . . .E{un). Recall that P^^'^ is defined by 
Using Lemma 6.5), we can express x'^\u(n) as ^ c{X)x^, where 

c(A) = c(Ai, . . . , Aat) = det[c(Afc - fc + j)]i<j,fe<Ar. 
Set fj{x}S) — c{xk + ])■ Since x'^ = x'^ ' dim^v A, with 

dimArA= n ^^^^^^^ = const • [] ((A^ - *) - (A, - j)) , 

we get the additional Vandermonde determinant det[(Afc — ky^^] — det[a;ir^]. 

□ 

Remark. Observe that the argument above and (jsj) imply that pT^-f {Cr) = 

det[fj{xk)\i<j,k<N- 

To state the next result we need slightly different notation. Let 

7''^^'''" ({4"^ : 1 < 71 < A^, 1 < fc < n}) = P''^^'^"(C^) 

if there exists a path r = (A'-^-' ^ . . . ^ A*^"-*) such that 

A(") = (4") + l,x(")+2,...,4")+n), 

and V^^'i {{x^k^}) = otherwise. 

Proposition 3.5. Let {a;^"-* : 1 < rt < A^, 1 < fc < n} &e arbitrary integers 
satisfying a;^"'' > x''^^^ for all n, k. Then 
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-pi+.i ({^.C^") : 1 < „ < AT^ 1 < fc < n}) = 

N-l 

const - J]^ det[(/)„(x|"\xj""^^^)]i<jj<„+idet[/j(x^-^^ 

n=l 

where x^^^i are virtual variable^^ and ((>„ is defined by 

{1 if x<y, 
1 X virtual^ 
otherwise. 



\l<i,j<N 



Proof. By the remark after 3.4 it suffices to prove that H^et [(/)„] acts as a 
indicator function, ft takes the value of 1 if A*^^^ -< A'-^' ^ . . . ^ A^^-' and 
otherwise, where 



If Ad) -< A(2) ^ . . . ^ A(^), so that > x^"^ > > a;^"^ > . . . > 

xl"'' > a;!'^^^'' for each n, then 

^ 10 it I < J. 

So det[0„] = 1 for each n. 

Conversely, suppose that H^et [(/)„] = 1, so that det ^ for each n. 
Notice that the matrix [(/i„(a;-"\ x^"'*'^')] consists entirely of zeroes and ones. 
Also notice that the number of ones in the fcth column is greater than or equal 
to the number of ones in the jth column for k < j. Additionally, if the («, j) 
entry is zero then so is the {i — l,j) entry. Since the determinant is nonzero, 
this means that no two columns are equal, so each column must have a different 
number of ones, so the (i, j) entry is 1 if j < i and if z < j. This says exactly 
that A(i) -< A(2) -< . . . ^ X^''^\ and each determinant in the product is equal to 
1. □ 



We can now prove Theorem 3.2 



Proof of Theorem 3.2 For computational purposes, it is actually easier to 
consider 

r if X < 2/, 

(f>n{x,y) := < 9~y X virtual, 
I otherwise. 

with mutually distinct 0„'s and then take 6'„ ^ 1. It is also convenient to 
denote 6*0 = 1. We will assume that \9n\ > \On-i\ > 1 for all n. Notice that 
det[/i(Xj^-')] only depends on the hnear span of /i, . . . , /at (up to a constant), 
so redefine 



^Onc can think of virtual variables as being equal to negative infinity. 
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fj{x) = - — : (p E{u)u ^pj^i{u ^)du, where 

^'"'^ J \u\=const 

N-1 

Pj.i{x)=^{9o-x)...{9;^x)...{eN-i-x)= l[ {9k ~x). 

k=Q,k^j-l 

The rest of the proof is a direct apphcation of Lemma 3.4 of [7J, where we 
use the notation — fj- 

Taking the Fourier Transform of we obtain 

<l>n{x,y) ^ ^ (f Fn{z)z''-y-^dz 

^(«l,n.)(^^y)^ 1 i Fr,,{z)...Fr,,-l{z)z^-y~Hz (8) 



where F„(z) = (1 — 6^^z)~^ and ni < n2- We also agree that = q if 

ni > 712. lu case a; is a virtual variable (which is denoted by virt), then 



9~y f dz 
■^|e„ii<|z|=const<|e„j+i| 2; - 

^KiiFni+i{9m) ■ ■ ■Fn2-^i{9m) 



This allows us to calculate the matrix M (cf. T, Lemma 3.4). In the 
following equation, r(ri, denotes the boundary of an annulus of radii ri < ri 
in the complex plane. 

^ 27rz/|„|=i 

^TTl Jr(r,l),r<|e,-i|-i i - Mt/.j_i 

=F,(0,_i) . ..FNM9,_i)Ei9r-_\)p^_^{9,.i)9,.i. 
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Notice that M is diagonal because = unless i = j. We have one 

more preliminary calcuation (cf. [7|, formula (3.22)): 

yez 

= (w-) { Fn{z)...FN-i{z)z''-^dz (f E{u)u-^pj-i{u-^)y{zu)-ydu 
= (^) i Fn{z)...FN-i{z)z-^dz(j> E{u)u-^-'^p,^i{u-')- 

V^TTiy Ji^i^i J\u\ = R>l 1-[ZU) 

= 7^ i Fniu-')...FN-liu-^)Eiu)u-^'''pj.i{u-')du 

27ri 

We can now calculate K according to f7^, formula (3.26). For ni < n2, 

= J2[M-%kKl-ki^l)^^''-'•"'Hv^rt,X2) 
fc=i 

= (f F,,^{u-^)...Fn_^{u-^)E{u)u-^-^' 

271-1 J\u\ = l 

P. O^^lFk{0k-i) ■ ■ . Fn,-ii0k-i)Pk-iiu-^) 
X > du 

^ Ffc(0fc_i) . . . FN^i{ek~i)E{e^l^)pk-ii9k-i)ek-i 

= ^ i Fr,,{u-^)...FN-i{u-^)E{u)u-^-^' 

xV O.-iP^-^ju ) 



k 

1 



Fn2{0k-i) ■ ■ ■ FM-i{0k-i)E{9f^_^)pk-i{0k-i)dk-i 



U 

27ri .l\u\=i 



X 



^^^^^ 6'„2 . . . eN-iE{ei_^) nr=o,i#fc-i(^' ^ ^fe-i) 

"2 a a "rr"i~i/ 



2:2 



fc=i+„i ~ ""^)^'n2 • • • (^N-iE{9k-i) nr=o,iy4fc-i(^' ^ ^fc-i) ^ 

and for ni > 712 the last sum is omitted. 

We can write the expression in parantheses as a contour integral that goes 
around all the 9j, so we get 



27^^; 7|„|=.-i>i7|,_i|=, (60- z)...{en,^i- z)E{z-^)z^+^^ (y-i - z) 
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assuming that |0„ — 1| < e for all n. Substituting u u ^ gives 

M V / / (gp ~ . . ■ (0^-1 - u)E{u-^)u^^ . ■ . ^^^^ 

There is also the term — (/)'^"i'"2^(a;i, 2:2) from ([s]), which equals 

~ V^iJ f\z\=const<l (1 - 9^1 Z) ... (1 - 0,-l_iZ) ^ 

if rii < and if ni > ri2- Finally, taking all the 9j to be 1 yields the result. 

4 Limits 

4.1 Limit Shape 

Represent A £ GT^r as a pair of Young diagrams (A+, A^). Figure [2] gives an 
example with A = (4, 2, 0, 0, -1, -3), A+ = (4,2),A- = (3,1). We have the 

Figure 2: The double lines show the boundary. 



following conjecture: 

Regard A S GT^r as random objects on the probability space (GT^r, ''^ ). 
As ^ 00, the boundaries of the two Young diagrams, scaled by N~^/'^ , tend 
to (nonrandom) limit curves. Both limit curves coincide with the limit curve 
arising from the Plancherel measure on symmetric groups. 

Our results strongly suggest that this statement holds, see §3.2. 

The conditions — (3^ — tell us that for fixed 7^ every row and column 
length grows sublinearly in N (see the end of §1). Furthermore, since 7* corre- 
spond to the area of the Young diagrams A^ (see §1), this suggests a scaling of 
iV-1/2. See Figured 

Furthermore, we see from Figure [l] that vertical segments of the boundary 
correspond to points in the configuration, while horizontal segments correspond 
to points not in the configuration. This implies that the first correlation function 
pi {x) (also known as the density function) corresponds to the density of vertical 
segments in the boundary. For example, in between the two curves in Figure [3j 
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Figure 3: A visual representation 
2^/7+ Vjv 




the vertical segments are densely packed, so pi{x) should converge to 1. Above 
the top curve (the boundary of A"*") and below the bottom curve (the boundary 
of A^), the horizontal segments are densely packed, so pi{x) should converge to 
0. We will see that this is indeed the case. 

Notice that near the edges of the Young diagrams (the boxes in Figure [s]), 
the probability of finding a vertical segment tends to or 1. This means that 
the vertical segments (or horizontal segments) become so rare that they occur 
infinitely far away from each other. In other words, for any fixed fc, the differ- 
ences A^ — ^k+i ^-'^'^ i-^^Yk ^ ('^*)fe+i both go to infinity as iV ^ 00. In fact, we 
find that A^ - A^^^ and (A±)'j. - (X'^^Yk+i are of order N^/^. The hmiting distri- 
bution of A^ — A^^^ or (A^)'j, — (A^)'j,_|_;^ normalized by N^^^ is referred to as the 
edge scaling limit. We will later prove that the well known Airy determinantal 
point process appears in the edge limits. On the other hand, if we zoom in at 
any other point on the limit curves, the behavior there is different. At these 
points, the differences between consecutive rows and columns stay finite. Their 
limiting distributions are described by the bulk limit. We prove that it coincides 
with the discrete sine determinantal process. The limit density function in the 
bulk predicts the limit shape. 

We should also consider what happens to the more general object - the 
corresponding measure on the set r of paths in GT (see §1). Consider two 
signatures on such a path at levels rii and n2 . If ni — n2 stays bounded then the 
bulk and the edge limits of these two signatures are indistinguishable (the local 
point configurations are essentially the same). However, as n2 — ni grows, we 
may see nontrivial joint distributions. It turns out that the proper level scaling 
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in the bulk is ni — n2 ^ \/N while at the edge it is ni — n2 ^ iV^/^. We will 
compute the corresponding scaling limits of the correlation functions later. 

It is also interesting to consider the case when the parameters 7* depend 
on A^. If 7* depend on TV in such a way that 7^7V ^ a > 0, then the areas of 
the Young diagrams stay finite. More precisely, we obtain two independent 
copies of the Poissonized Plancherel measure for symmetric groups. 

Additionally, consider what happens when 7^ depend on TV in such a way 
that 7+/iV ^ a > and 7~/iV^6>0as7V^cxD. The Young diagrams are 
now scaled by 7V~^ The new hypothetical limit shape depends on the values 
of a and b. See Figure |4] 

The edges of the limit curves correspond to the real roots of a fourth degree 
polynomial 

1\ / 1\' / J \ ' 



Qa,6(z) =P0+Pl 2j +P2 (^^+ 2j +?'3(^2+-j +16^,. , 2 

po = 1 - 12(a + 6) + \{o? + l?) + 184a& - 256a6(a + 6) + 64a6(a - 6)^ 
= 8(6 - a)(7 - 2a - 26 + 16a6), 
P2 = 8(2(a + hf - 10(a + 6) - 1), = 32(6 - a). 
The expression Qa.6(c) is the discriminant of a simpler polynomial 

Ra,hA^) = + (6 - c - l)z^ + (c + a)z - a. 

For small a and 6, Qa.b has four real roots. As a and 6 increase, two of the real 
roots become closer until they merge into a double root. For larger values of a 
and 6, Qa.b{z) has two real roots. 

We will be able to find what values of a and 6 lead to Qa.b having exactly 
three distinct real roots (the middle root is a double root). This corresponds to 
the situation when the two limit curves just barely merge (see the middle image 
in Figure |4|. The correct scaling there is to let (A*)J — {X^Y^^i '-^ iV^/** and 
ni — ~ N^^'^, which results in the Pearcey determinantal process appearing 
in the limit. At the other edges, letting - A^+i iV^/^ or (A='=) ■ - (X^Yi+i 
and ni — 712 N'^l^ results in the Airy process appearing. Away from the edges 
we still observe the bulk limit. 

We now proceed to computing the (scaling) limits of our determinantal point 
process corresponding to the limit regimes described above. 

4.2 Limits with 7=^ oc 

Introduce the kernel J on M+ x Z by 



{s,x;t,y) 



1 V r r e" '-t^-^ '+^^" dudw 



2711 ITT w — u w^'^^u y 



where the w and u contours go counterclockwise around in such a way that 
the w-contour contains the u-contour if s > i, and the w-contour is contained 
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Figure 4: The limit curves for various values of a and h. The top curve occurs 
when a = = , the middle curve occurs when a = b = | , the bottom 
curve occurs when a= \,b= | . 
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in the w-contour if s < i. This kernel for s = t is equivalent to the discrete 
Bessel kernel KbcssoI, which appears when analyzing Plancherel measures for 
symmetric groups (see e.g., §2.4 of [25). Additionally, J is a special case of the 
kernel (|TT], (3.3)) corresponding to 9{t) = e~^*. 

Theorem 4.1. Let xi, . . . ,Xk be finite and constant. Let ni, . . . ,nk and 7^ 

depend on N in such a way that rij/N tj and •y^N a > 0. Then as 
N ^ 00, 

det[K{ni, x^; nj,Xj)]i<ij<k det[:S{ati,Xi;atj,Xj)]i<i^j<k 



Proof. We use the integral representation for the kernel in Theorem 3.2 

We first focus our attention on the double integral in u and w. Since the 
integrand is holomorphic everywhere except ai u ^ 0, w ~ I, w =^ u and w = 0, 
we can deform the contours of integration as shown in Figure [5] 



Figure 5: Deformation of the contours. 




As ^ 00, the integrand converges to for \w\ large enough because 
|1— ti;| 3> |1 — u|. Therefore we can ignore the outer half of the w contour. 
Then the contours of integration can be deformed to |m| = a/N and \w\ = 2a/N. 
Making the substitutions u' = Nu/a and w' = Nw/a, the double integral is 



1 / / eT^"' u'a/N u'^. (1 _ y'a/7V)"* du' dw' (N 



(27rz)2 7|„,|^i J\w'\=2 ei^^'~^^/''+^'^'°-/^ u;'^j+i(l - w'a/NYo w' - u' \a 

1 f f ^u'-^-aUu'-w'-^+at,w'+Oil/N) ^^^^ , _^ 

(27rz)^ /|„,|^i/|„,|^2 w'-u' u' ^-w'^^+^ 
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When taking the determinant, the term {N/y/a)^^~^' cancels. This gives the 
result. 

Remark. Comparing this result to [8j, we see that the distribution of A+ 
converges to the Poissonized Plancherel measure for the symmetric groups. By 
the symmetry (A^ <-> A^, 7^ <-> 7^) the same is true for A~. On the other hand, 
a similar contour integral argument to the above shows that if (n^, xf, rij, —rij — 
Xj — V) — > as A'^ ^ 00, which implies that A"'' and A~ are asymptotically 
independent. □ 

4.3 Bulk Limits with 7=*" fixed 

To state the next result, we need a definition. Given a complex number in 
the upper half plane, define 

S,^{U - tj-x, - X,) ^^■J_^ u^^-^'^-^e^'^-'^^du. 

liti > tj, then the integration contour crosses (0, 00) but does not cross {—00, 0). 
If ti < tj, then the integration contour crosses (— oo,0) but not (0,cxd). This 
kernel is one of the extensions of the discrete sine kernel constricted to |5] . A 
similar kernel appeared in . It can be seen as a degeneration of the incomplete 
beta kernel, see Section [441 

The main theorem of this section is the following: 

Theorem 4.2. Let xi, . . . ,Xk and ni, . . . ,nk all depend on N in such a way 
that Xi — Xj is constant, {uj — N)/y/N —> tj £ R and Xi/VN c £ R for all 
l<ij<k. Write z+ for (c+ ^Jc^ - 47+)/2. Then 

lim det[K{n^,x^;nj,xA]i^i.<k = 

N—>-oo 

0, c > 2^7+, 

1, c < -2v/7+, 
det[Sz+{tt - tj;xi - a;j)]i<i j<fe, -'2\/r^ < c< 2^/j+. 

Remark. Theorem |4.2| only makes a statement about the behavior around 
the top limit curve in Figure [sj If we replace Xi with —Xi — — 1 and 7+ with 
7~, then by symmetry the same statement holds for the asymptotics around 
the lower Young diagram. 

Corollary 4.3. Let pi{N,x) be the density function of V]^ . Then lim pi{N,aN+ 

N^oo 

f3N^/^) equals 

• 0, if a > or a < ~1 or a = 0,/3 > 2^7+ or a ^ -1,/? < -2^7+. 

• i < a < or a = 0,/3 < -2^7+ or a ^ -1,(3 > 2y/j+. 



2V7 



if a^O, -2^7+ < /? < 2^7+ or a = -1, -2^7+ < 



f3 < 2y/^ 
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Proof. The arguments are similar to those used for the analysis of Plancherel 
measures for the symmetric groups in |20j . 

For reasons that will later become clear, it is more convenient to analyze 

N ' ^ ' ^ K{ni,Xi;nj,Xj). When taking the determinant 



det[-\/7V ' ' ' ' K{ni,Xi;nj,Xj)], 

the term \/N ' ^ ^ ' cancels out. 

We use the integral representation for the kernel in Theorem |3.2[ The con- 
ditions Hi > Hj and Ui < rij translate to ti > tj and ti < tj, respectively. 

Just as in Theorem |4.1| we can deform the contours of integration as shown 
in Figure |5] 

As — > cx), the integrand converges to for |ti;| large enough because 
|1 — lyj 3> |1 — w|. Therefore we can ignore the outer half of the w contour. Then 
the contours of integration can be deformed to |u| = 1/VN and |it;| = 2/^/N. 
Making the substitutions u' = ^/Nu and w' — ^/Nw, the double integral is now 



(27rz)2 ij„|^i/y]v f\w\=2/VN eT+--'+^-- w^^+^e^^ '"(i—) w-u 

7V7+ 



' (27rz)2 l\w'\=1 e7+«''-'v7V+7-u)'/v7V yjlXj + \^n, ln(l-i«'/v7V) w' -u 

If f g7]v(7+«'"'+ciogu'-u'+o(i/v7v)) du'dw' 
' (27ri)2 J\w'\=2 e^i-r+w'-^+ciogw'-w'+oii/VN)) w'{w' - u') 



In general |e^| = e^^ , so consider the real part of the function in the expo- 
nent, A{z) = 7+z~i -I- clog z - z. Note that A'{z) = at z± = f ± ^LfE^jl^ 

The basic idea of the rest of the proof can be summarized as follows. The 
term 1/7+ ^ ' creates a e^^'^~'^'°s term in both the numerator and de- 
nominator. So it is equivalent to analyze 5R(A(z) — clog \fy^) = 3?(A(z) — 
A{z+y). We deform the u and w contours in such a way that 3fi(A(u)— A(z+)) < 
and ^{A{'w) — A{z+)) > 0, which will cause the integrand to converge to as 
N ^ Q. However, the deformation of the contours causes the integral to pick 
up residues at u = w. These residues occur on a circular arc from z„ to z^. If 
c = 2a/7+, then z+ = z- > 0, so the arc consists of a single point. As c de- 
creases, moves counterclockwise around the circle \z\ = \px^ while Z- moves 
clockwise. This means that the arc becomes increasingly large as c decreases 
from 2-^/7+ to — 2-^/7+. When c = —2-^7+, then z+ = z_ < 0, so the arc has 
becomes the whole circle around the origin. 

We then need to consider 

1 / 

<P -, ^ 

2« /|.|=.<i (1 - ^)"^-"' 
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which occurs when Ui < rij. The expression for the residue at u — w has 
the same integrand. With the minus sign, the integration contour for z goes 
clockwise along a circle around the origin. Therefore it will cancel the circular 
arc from z„ to This explains why the integration contour in S^^ crosses 
(0,cxd) when ti > tj and (— oo,0) when ti < tj. 

Case 1: -2yfi+ < c < 2yfi+ . Observe that 3fJ(A(z) - A{z+)) = for all 
\A — l-^il — IVt^I- Also notice that A{z) — A{zj^) has a double zero at z+ and 
Z-. See Figure [6j 

Figure 6: On the left is 5ft(A(z) — clog Vt^), where the black regions indicate 
5J < and the white regions indicate 3? > 0. 




If the contours of integration are deformed as shown in Figure [6] then 

„V]V(7+ti'~^+clogti'-M'+0(l/\/]V)) 

^ ^ -^0 

g^s/N(^+w'-^+c\og,w'-w'+0(l/\fN)) 

as N oo. The integral thus approaches zero, except for the residues at u ^ w. 
So y/N ' ^ K{ni,Xi;nj,Xj) converges to 

If ti > tj, then the integration contour crosses (0,cxd). If ti < tj, then the 
contour crosses (— oo,0). 

Case 2: (? — 47+ > and c > 0. Deforming the contours of integration as 
shown in Figure |7] the integral becomes zero. The contours do not pass through 
each other, so no residues appear. So \fN ' ^ K{ni,Xi;nj,Xj) — > if > tj. 
This means that det[K(ni,Xi;nj,Xj)] — > 0. 

Case 3: (? — 47+ > and c < 0. Deform the contours as shown in Figure [s] 
Since the w and u contours pass through each other during the deformation, the 

integral picks up residues at u — w. So if > tj , then \/iV ' ^ K(ni , Xi ; Uj , Xj ) 
converges to 
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Figure 7: Again, the figure on the left shows ^{A{z)—A{z+)), with black regions 
indicating 5i < and white regions indicating SR > 0. 




2m 



w 



2m 

(Xj — Xi)\ 



-ie(t.-t.)«'du; 



If ti < tj , then there is the integral in z, which cancels with the residues at 

u = w, so Vn ' ^ K{ni,Xi;nj,Xj) converges to 0. This means that the matrix 
[K{ni,Xi; nj,Xj)] asymptotically has ones on the diagonal and zeroes below. So 
det[K{ni, Xi; rij , Xj)] converges to 1. 

□ 

Remark. It is natural to ask what happens when Xi/y/nl do not all converge 
to the same real number. When this occurs, the determinant dct[K{ni, Xi; nj,Xj)] 
factors into blocks corresponding to distinct values of lim .x^/y^. Probabilisti- 
cally, this means that the probability of finding a vertical edge becomes inde- 
pendent in different parts of the boundary. 

4.4 Bulk Limits with 7^ oc 

We now let 7^ depend on N in such a way that j^/N — > a > and 7~/iV 
b > as N ^ 00. Before we can state the result, some preliminary definitions 

and lemmas are needed. 

For a,b > and c G M. recall that 

Ra,b,c{^) = -i>z^ + (6 - c - l)z^ + {c + a)z - a. 
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Figure 8: On the left is 3?(A(z) — clog \/7+), where the black regions indicate 
5i < and the white regions indicate Ji > 0. 
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Lemma 4.4. (1) The cubic polynomial Ra,b,c{z) has a multiple root iff c is a 
root of Qa^b{z), where Qa.h{z) is defined in ^3.1. 

(2) Let qi < ■ ■ ■ < Qm be the real roots ofQa^- If qi < c < q2 or qm-i < c < 
qm then Ra,b,c{z) has a pair of complex conjugate roots. 

Proof. (1) In general, a polynomial has a multiple root iff its discriminant is 
zero. The discriminant of Ra,b.c is exactly (5a,b(c)/16. 

(2) A cubic polynomial has nonreal roots iff its discriminant is negative. 
Since Qa,b diverges to +oo in both directions, Qa,b{z) is negative for qi < z < q2 
and qm-i < z < qm- □ 

Lemma 4.5. The polynomial Qa,b has a double root at cq iff a, b and cq satisfy 
the equations 



izo-iy 



1 

Jzo^' 



Co 



{zo-ir 



(9) 



for some Zq G 



Proof. Since Qa,b{co) is the discriminant of Ra,b,coj Qa,b{z) has a double root 
at Co iff Ra,b,co{z) has a triple root. For any zq, R has a triple root at zq iff 
R{zo) = R'{zq) = R"{zo) = 0. This gives three linear equations in the three 

variables a, b, and cq, which can be solved explicitly. □ 

Remark. Wc have a, 6 > iff zq < 0. Then — 1 < cq < 0. 
One more definition is needed before we can state the main result of this 
section. Let B be the incomplete beta kernel defined by 
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Figure 9: This figure shows the equations in with a plotted on the horizontal 
axis and b plotted on the vertical with parameter zq. 




where the path of integration crosses (0, 1) for A: > and (— oo, 0) for A: < 0. The 
incomplete beta kernel has been introduced in [23]. It is one of the extensions 
of the discrete sine kernel of [5^ . 

Theorem 4.6. Let /N a and "/^ /N b for positive real numbers a and 
b. Also let Xi, . . . ,Xk and ni, . . . , n^. depend on N in such a way that rij — rij 
and are Xi — Xj constant, Uj/N — s- 1 and Xj/N — s- c for all I < i, j < k. 
Let qi < ■ ■ ■ < Qm denote the distinct real roots of Qa,bix) (m can be 2,3, or 
4). Additionally, assume Qa.b{c) ^ 0. Let be a root of Ra.b.c{x) such that 
> (cf. Lemma 4-4)- ^/'ti = 4, then 



det[K{ni,Xi; nj,Xj)]i<i^j<k 

If m = 2 or 3, then 
det[K{ni,Xi;nj,Xj)]i<i^.j<k 



0, 

det[B(nj 
1, 

det[B(ni 
0, 



c < qi, 

nj,Xj - Xi; z+)]i<i,j<fc, gi < c < q2, 

q2<c< q3, 

nj,Xj - Xi; z+)]i<i,j<k, qs < c < q^, 

c > qi. 




Tij , Xj Xi , -^-j- )] l<i,j</c ; 



gi < c < qra, 
c > qm, 



Proof. The double integral in the correlation kernel of Theorem 3.2 asymptoti- 
cally becomes 



1 

27ri 



gAf(an-i+6«+clog(u)+log(l-«)+0(l/Af)) dudw 
jAr(atu-i+6u.+clog(u)+log(l-ji)+0(l/Af)) y;^^ — w) 



where the contours are over \u\ = r and |w — 1| = e < 1 — r. So we can 
perform a similar analysis as in Theorem |4.2[ except with a more complicated 
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Figure 10: The shaded regions show c) — c)) < 0, while the white 

regions show > 0. The first row corresponds to c < qi, the second row 
corresponds to gi < c < q2, the third corresponds to ^2 < c < (73, the fourth 
corresponds to 93 < c < 94, and the fifth corresponds to c > 54. 
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A{z) = az~^ + bz + clog(z) + log(l — z). For this proof, it is actually more 
convenient to write A(z; c) in place of A(z). 

First we find which values of c correspond to the edges of the hypothetical 
limit shape in Figure [4] These are the values of c such that A{z; c) — A{zo; c) 
has a triple zero for some zq G C. Requiring A(z; c) — A(zq; c) to have a triple 
zero at z — zq is equivalent to requiring A'{z;c) to have a double zero at 
z = Zq. Multiplying the equation A'{z;c) = by z^(l — z) gives the equation 
Ra,b,c{z) = (Note that i?a.&,c(0) = —a and Ra,b.c{^) = —1, which are both 



nonzero). By Lemma 4.4 Ra.b,c has a double zero iff c = gi, . . . , Qm- 

Now we need to determine how to appropriately deform the contours. The 
analysis here is almost identical to that of Thereom |4.2[ We want to find nonreal 
values of zq such that A{z; c)—A{zo; c) has a double zero. This reduces to looking 
for nonreal roots of Ra,b,c{z) in the upper half-plane, which we have defined to be 
z+. As can be seen from Figure |4] there are potentially five different regions of 
behavior for the bulk limits. The corresponding behavior of ^{A{z\ c) — A{zj^ \ c)) 



is shown in Figure 10 (These are computer generated figures for specific values 
of parameters, however, it is not hard to prove that similar figures arise for any 
values of the parameters in the corresponding domains. An example of such an 



argument can be found in the beginning of the proof of Theorem 4.7 below.) 
The arguments of Theorem |4.2| are again applicable here, except with the new 
definition of □ 

4.5 The Pearcey Kernel as an Edge Limit 

We now find the edge limit at the point where the two limit curves in the 
middle figure in Figure[4]just barely merge. In this case, we analyze the limiting 



behavior of K/^ from Corollary 3.3 instead of K, which corresponds to the fact 
that we consider the limit of the point process formed by columns of A* rather 
than by their rows, see Figure [T] 

Theorem 4.7. Fix Zq < and let a, b and Cq satisfy equations Let 7+/iV — > 
a and 7"/-/V — s- 6 as N ^ oo. Let ni, . . . , depend on N in such a way that 
{uj - N)/^/N 2tj eR as N ^ oo. Set C ^ (zq - l)\zo\-^/'^ < 0. Define 

i -^^t 
1-zo'' 

and let xi, . . . ,Xk depend on N in such a way that 

^175 '^^^^ 

as N oo. Then as N oo, 

det[-C~'^N'^^'^KA{ni,Xi;nj,Xj)]i<i^j<k -> det[P{ti, Si;t.j, Sj)]i<^,j<k 
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wht 



1 

27ri 



1 

27ri 



^w'^ -u^ +tiU^ -tjW^ + SiU-~SiW dudw 

u — w 

2 ^ 



^2n\U-tj\ 



exp 



jsj - g») 

dudw 



u ~ w 



U > tj 



— ^7 



(10) 



where u is integrated from — ioo to ioo and w is in tegr ated on the rays from 
ibooe*'^/^ to and from ±ooe~"/* to as in Figures\ll\ and 12 



Figure 11: The contour for u. 



Figure 12: The contour for w. 



The kernel P{ti, Si, tj, Sj) is called the Pearcey kernel and it was previously 
obtained in [J , US] , [II] , [M] , [H] . 

Proof. The argument is similar to the proofs of Theorems |4.2| and |4.6| It is 
convenient to let A{z; c; d) denote az^^ + &z + clog z + d\og{\ — z). Then the 
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double integral in the correlation kernel of Corollary 3.3 becomes asymptotically 



1 
1 

27ri 



^N{au-^+bu+(xi/N) log ti+(ni/JV) log(l-M)+0(l/A')) 



dudw 



gN{aw-^+bw + {xj/N)logw+{7ij/N)\og{l-w)+0{l/N)) uj{u — w) 

(11) 

^N{A{u;x,/N;n,/N)+Oil/N)) ^^^^ 



^N{A{w:x,/N;n,/N)+0{l/N)) _ 

Multiplying the integrand by the conjugating factor 

,~aNzZ 



-bNzo 



e 

Q-bNzo ' 



(12) 



^-NA{zo;x,/N;n,/N) ^ ^-x, (^l _ zo)~' 
g-NA{zo;Xj/N,nj/N) ~ ^-Xj ^-^ _ Zq)""^' e~°-^^ 

which cancels when taking the determinant for correlation functions, allows us to 
consider A{z; Xm/N; nm/N)—A{zQ\ Xm/N; Um/N) instead of A(z; Xm/N; Um/N). 



Figure 13: The figure on the left shows ^{A{z; cq; 1) — A{zo; cq; 1)), with black 
regions indiciating 3? < and white regions indicating 3? > 0. 





Deform the contours as shown in Figure 13 Let us show that these contours 
exist. We know that the level lines only intersect at zq (the only critical point 
of the function A{z] Cq; 1) — A{zo; cq; 1), since A'{z) = —b{z — ZQ)^z~'^{l — z)~^), 
and they are symmetric with respect to the real axis. Restrict ''St{A{z; cq; 1) — 
A(zo;co;l)) to the real axis. For = e small, the main contribution to 
Si(A(x;co; 1)) comes from the term ax~^ . So ^{A{x;co; 1)) is positive at x = 
e > and negative at a; = e < 0, so the level lines cross the real axis at 0. For 
X ~ 1 — e with e small, the main contribution to ^{A) comes from the term 
log |1 — x|. This implies that 3?(A) is negative a: = 1 — e, so the level lines cross 
the real axis somewhere between and 1. For large x, the main contribution to 
5R(A) comes from bx, so ^{A) is positive for large x. Therefore the level lines 
cross the real axis at a third point. Since A'{z) — ~h{z — zqYz~'^{\ — z)~^ is 
positive for z < zq, negative for z G {zq, 0) U (0, 1), and positive for z > 1, the 
levels lines can not intersect the real axis at any other point. 
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For a fixed x <C 0, the main contribution to ^{A{x)) comes from bx, so 
^{A{x)) is negative. However, as y increases, ^{A{x+iy)) goes to +00, since the 
main contributions come from colog jx + ij/j + log \ l — x — iy\, and cq > —1. This 
means there must be level lines going off to infinity. Restricting ^{A{z; Cq; 1)) to 
a circle \z\ = R ^ 1 shows that these are the only level lines that go to infinity. 
Indeed, note that ^{A{z; cg; 1)) > if z = R, and as z moves counterclockwise 
around the circle, the main contribution to the changes in ^{A(z)) comes from 
bz. Thus ^{A{z)) decreases as z moves counterclockwise around the circle in 
the upper half-plane, so the circle can intersect at most one level line in the 
upper half-plane. 

In the upper half-plane, there are four level lines coming from the critical 
point zq. We know that three of these lines cross the real axis, while one of 
them goes off to infinity. Since they can only intersect at zq, the only possibility 
is a picture as shown in Figure [13] This justifies the existence of the contours. 

These deformations cause the kernel to pick up residues at u = w. The 
expression for these residues is 



1 



-dz (13) 



27ri J (1 - z)-" 

where the integral goes around a circle |z| < 1. If < rij, then expression 
(13) cancels with the z-contour in expression ^. If rii > Hj, then explicitly 



evaluating the integral yields 



X 'i Xi 



The binomial can be approximated by the deMoivre-Laplace Theorem. For 
large N, 

\Xj — Xi) 



1 / {Sj - Sj) 



, : exp , „ , 



^2 



So when ti > i^, we obtain the extra exponential term in equation (10 1. 

For large values of N , all the contributions to the double integral come from 
near the point zq. Taking the Taylor expansion around zq yields 

(^^;^o+^ + ^;l + ^j-A(^zo;co + ^ + ^;l + ^ 

= s^z' + tMf-{^'Y + o{i) 

where z' — z'f^'^Q^'^N^/^{z — zq). This suggests the substitutions 

u' = z^\-^N^/^{u - zo), w' = Zq\-^N^'^{w - zo). 
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By making these substitutions, we are zooming in at the point zq in Figure 13 



Then u' is integrated as shown in Figure 11 while w' is integrated as shown in 
Figure [121 



The exponential terms in expression ( 12 1 converge to the exponential terms 

For large N, the contribu- 



in dlOl). The term turns into zoCiV~'*^¥^ 



tions to the correlation kernel become focused around zq , so the extra w in the 
Zq ^ . The proof of Theorem 



denominator becomes 



4.7 



is complete. 



□ 



4.6 The Airy Kernel as an Edge Limit 

Before stating the main result, some definitions are needed. 
Let Ai{x) denote the Airy function: 



Ai{x) 



1 

2^ 



This integral o nly c onverges conditionally. Shift the contour of integration as 



shown in Figure 14 Along this contour, the function e** is real and decreases 



superexponentially. 



Figure 14: A better contour for the Airy function. The contour goes from 

QQg57rj/6 to e""'/^. 




Define the extended Airy kernel A to be 

_ //o°° e-^(-i--^)Ai((7i + \)Ai{a2 + X)dX if n > T2, 



A{Ti,ai;T2,a2) 



- e-^(''i"^2Uz(tTi + X)Ai{a2 + X)dX if n < ra. 

(14) 



It was first obtained in in the context of the polynuclear growth model. 
There is a useful representation for ^ as a double integral. 

Proposition 4.8. ( fT5f . ^2.2) Let vi, V2 satisfy V1 + V2 + T1—T2 > 0. //ti > T2, 
then 

/ \ \^ f f gicriM+lcr2to+i(tu^+«^)/3 

A{Ti,ai;T2,a2) ^ [—.] / / —, ■ — :—dudw. 
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// Ti < T2; then 

AiTi,(Ti;T2,(T2) = [ —] / / — — —dudw 



v/MT2-n) "V 4(r2-ri) 2^'^ -m^--/'i2 
The double integral from Proposition |4.8| can be rewritten as 



exp (tio-1 - r2cr2 - + -t| - (cti - )u + ((72 - t|)w 



-TlU^+T2W^+ -(U -W-^)) . (15) 

6 U — W 

Indeed, just as we deformed the contours of integration for Ai{x), we can deform 
the contours of integration in Proposition |4.8[ The u-contour can be taken over 
ivi + ooe^^^^^ to ivi to ivi + e^^l^ , while the tu-contour can be taken from 
iv2 + ooe^'^'/^ to W2 to iv2 + e^^l^ . Integrating along these contours also allows 
for the possibility of i^i + 1^2 + — T2 = 0. If we further make the substitutions 
w = —ivJ + i^2* and u — iv! + i>ii, then the double integral becomes 

J y exp ( - i^iCTi - 1^2(^2 + l^l^l + ^i^l - (cti - I^f)u + (cr2 - 1^^)11! 

9 , 1 , c, dudw 

+ ViU-^ + Z^2W + -{U-^ - W'^)) . 

3 -T2 + Tl + Vl + V2 + U - W 

where u is integrated from cx)e~'^*/'^ to to cx)e'^'/^ and w is integrated from 
gQg47ri/3 Q g^^g27ri/3^ Taking vi = — Ti and j/2 = T2 turns the double integral 
into (15 1. Writing the double integral in this form is useful when proving the 
following result. 

In the next statement, let Qa b be the same polynomial as in §3.1, see also 
§3.4. 

Theorem 4.9. Let "/^/N — > a,^^ /N b for positive real numbers a and b. 
Let ci be a root ofQa.b{z) and zi be the double zero of Ra,b,ci{z)- Let ni, . . . 
depend on N in such a way that 

-j^ ^ e M as TV ^ 00. 

Let tj = tjZi{l — zi)~^ and let xi, . . . ,Xk depend on N in such a way that 

Xi - ciN - i,N^/^ 
, > Si eR as N ^ 00. 

If Ci > or ci < —1, set K, — K. Otherwise, set K, = A'a. Then as N — > 00, 
det[\zipy^\N^/'^K{ni, Xi, nj,Xj)]i<ij<k ~* det[^(Ti, ai; Tj,crj)]i<ij<k- 
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Here, ps denotes the constant 



3a ci 



and 

Remark. The statement may seem a bit cryptic. Let us explain it in words. 
There are (potentially) four edge points as seen in Figure |4j We consider K for 
the first point (when ci > 0) and the fourth point (when ci < —1), which means 
that we look at the largest rows of and A^. For the second and third points 
we consider K/^ , which means that we look at the largest columns of and A^ . 

1 /3 

For the second and fourth points, det[0iP3 JC] det[yl], while for the first and 

third points det[—Zip^^^]C\ —^ det[^]. At the second and fourth points 2:1^3^^ is 

1/3 

positive, while at the first and third points 21^3 is negative. This corresponds 
to the fact that in order to obtain the Airy process we need to flip the sign of 
particles at the lower edges of A+ and A~ (the second and fourth edge points, 
respectively) . 



Proof. This proof is similar to the proof of Theorem |4.7| so some of the details 
will be omitted. 

Once again, let A{z; c; d) denote az~^ + bz + c\ogz + d\og{l — z). Multiplying 
by the conjugating factor 



allows US to consider A{z; x,n/N; nm/N)—A{zi;Xm/N\ rim/N) instead of ^(z; Xm/N\ rim/N). 
The Taylor expansion yields 

N [Al z;ci + —-p; + — 1 + — - A zi;ci + — + — jt^t; 1 



7V1/3 7V2/3' iVVsy \^ i'-i ^1/3 ^2/3' JSll/3 

3 2(p3)^/-^(zi - lyzi (p3)V3zi 

where z' = {p-^) ^/^ N^/^{z — zi). The contours of integration for u and w axe 
shown in Figure [isl Now let u' = {p^Y/^N'^/^{u- zi) and w' = {psY/^N^/^iw- 
z\ ) . Just like in the proof of Theorem |4.7[ the Taylor series gives rise to the 
exponential terms injisj In addition, the term becomes 7V~^/'^P3 while 
the extra w in the denominator becomes We break down the following 

analysis into cases. 

Case 1: ci > 0. This corresponds to the fourth row in Figure [T5| and the top 
edge point of Figure [sj In this case, p^ is negative, so the contours for u' and w' 



agree with the contours in expression (151. Since < zi < 1, this implies that 



tj — > if tj — ti > 0. Since rij > rii translates to tj > ti, this means that 
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Figure 15: The left column shows lR.{A{z; cq; 1) — A{zo; cq; 1)), with shaded re- 
gions showing 3? < and white regions showing 3? > 0. The right column shows 
the local behavior around zi. The first row occurs when ci is the smallest real 
root of Qa.b, the second row when ci is the second smallest real root, and so 
forth. If Qa.b has only two real roots, the middle two rows do not occur. 




u w 



1 
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be assumed positive if nj > n^. Therefore the integral in z from 
expression Q can be written as 

Uj — fij + Xj — — 1\ / Uj — Hi + Xj — Xi\ Uj — rii 

"X J X> 1 J \^ J ^ / 3 ~ ' I X J X 

Using the Laplace-Demoivre Theorem shows that 



-iVl/3 



(1 - zi)"--«^ - ^^ 



:expf- (\7^|)'(p-^f ). (16) 



^27:\zi\itj-U) "V 2|zi| |i 



Taking the last term in Proposition 4.8 and multiplying by exp(— tiCTi + T2<T2 + 
- 5^1) yields 

We have seen that 

\zipI^^\N^/^ ^^_^^^^^_^^^ K{ni,Xi]nj,Xj) -> 

exp(-ricri + T2(J2 + - ^t|)^(tj, (j^; Tj,(Tj), (17) 



which gives the result. 

Case 2: ci < —1. This corresponds to the first row in Figure [T5| Here, 
zi > 1 and P3 > 0. Making the deformations gives residues at u = which 
can be written as 

~2^/I.-i|=.<i (1-z)".-"' 

If Hi > jij, then these residues are zero. If rii < Uj, then U < tj, which implies 
Xi > Xj, so the integral in z from expression Q is zero. So when rt,; < Uj, the 
extra term can be written as 



fii — rt, — 1 / \ni — n 



Using Laplace-Demoivre, this binomial converges to 16 So expression (17 1 
holds. 

Case 3: — 1 < ci < 0. If a and b are small enough, then Qa.b has two roots 
between —1 and 0. The second row in Figure 15 corresponds to the smaller 
root, while the third row corresponds to the larger root. In the second row 
is positive, while in the third row is negative. In both rows zi < 0. 
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Making the deformations gives residues at u = w, which can be written as 



1 



^Xi—X-i — l 



-dz. 



If rii < Uj, then this expression cancels with the z-integral in expression (jsj). If 



n,; > 71 j , then the extra term can be written as 



Once again, this converges to expression (16). So expression (17 1 holds. 



□ 
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